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ABSTRACT 

We present a set of 144 galactic chemical evolution models applied to a Milky Way analogue, 
computed using four sets of low+intermediate star nucleosynthetic yields, six massive star 
yield compilations, and six functional forms for the initial mass function. The integrated or 
true yields for each combination are derived. A comparison is made between a grid of mul¬ 
tiphase chemical evolution models computed with these yield combinations and empirical 
data drawn from the Milky Way’s disc, including the solar neighbourhood. By means of a 
methodology, applied to the results of these multiphase models, the best combination of 
stellar yields and initial mass function capable of reproducing these observations is identified. 

Key words: Galaxy; abundances - galaxies: abundances - stars: abundances - stars: mass- 
loss - stars: supernova 


1 INTRODUCTION 

Chemical evolution models (CEM) were developed initially to un¬ 
derstand observations such as the local metallicity distribution of 
G/K-dwarfs and the radial gradient of abundances through the disc 
of late-type spirals, including the Milky Way Galaxy (MWG). The 
basic framework for a CEM involves a volume of a galaxy within 
which gas is assumed to flow, both inwards via infall and radial 
flows, as well as outflows; an adopted star formation prescription, 
coupled with an initial mass function (IME), then allows the calcu¬ 
lation of the production rate of stars of a given mass, supernovae, 
and the ejection rate of nucleosynthetic products back to the in¬ 
terstellar medium (ISM). The latter is often characterised via the 
use of stellar yields and the in tegrated or true yields, concepts first 
introduced bv iTinsle^ lll980li . While the infall and star formation 
rates (SER) are essential to create and maintain a certain radial 
abundance gradient, the IMF and the stellar yields define the abso¬ 
lute level observed in a region. It is therefore critical to understand 
the origin of the elements and, in particular, the specific stars from 
which individual chemical elements originate, the quantity returned 
by each said star, and the timescale for their ejection back into the 
ISM. 

Since the seminal work of [Burbidge et al.l (1 19571) much has 
been done to improve our understanding of stellar nucleosynthesis. 
Early work focused on stars with metallicities similar to that of the 
Sun, with contemporary work now concerned with spanning the 
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full range of metallicities encountered in nature. However, code- 
to-code differences still result in substantive differen ces in the pre¬ 
dicted stellar yields across these mass ranges (e.g. iGibson et all 
Il997h . Due to the quite different end-lives of massive stars, as op¬ 
posed to low -t- intermediate-mass stars, stellar evolution codes have 
typically separated their applicability to either those which end 
their lives as Type II supernovae (SN-II) or those which end their 
lives as white dwarfs. 


For massive stars, the total yields of elements are usually 
provided for those originating from the supernova explosions or 
those originating from pre-explosion stellar winds; only rarely 
are both provided, self-consistently. The most frequently used 
set of massive stel l ar yie lds (hereinafter mas) has been that of 
IWooslev & Weav^ d 19951 hereafter WOW); to the elements pro¬ 
duced in SN-II (for metallicities spanning Z - 0 to Zff), W OW 
added the pre-supernova yields of IWooslev & Weaver! (Il986l) . but 
did n ot include the contribution from pre-SN-II stellar winds. 
Later. fPortinari. Chiosi. & BressanI d 19981 hereafter PCB) provided 
massive star yields for a range of metallicities, but now tak¬ 
ing into account both the pre- and post-explosion elemental re¬ 
turn rates, including the stellar winds and the subsequent effect 
of this mass loss on the evolution of the sta r and on the ejec¬ 
tion of the supernova explosion. Mor e recently. iLimongi & Chieffa 
d2003h and Chieffi & Limongi ( 2004 hereafter both sets referred to 
as CLI), ILimongi & Chieffil d2006b and [Limongi^^Qiied 1 2012, 
hereafter both sets referred to as LIM), Kobavashi et al.l 1 200C, 
her eafter KOB ) and Rauscher et al.l d2002ih Frbhlich et all d200d) . 
and iHeger & WooslevT i 20101 her^fter the joined sets referred to 
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as HEG) have calculated new massive stars yieldfl We will use all 
these sets in this work. 

For the low and intermediate mass st ar (hereinafter 
Jim) yields, besides the seminal work of iRenzini & Volil 
(Il98lh . where the effects of convective dredge-up and the 
so-called Hot Bottom Burning processes were taken into 
account, more recent yield compilati o ns have been pro- 


Forestini & Charbonnell h 99' 


Ventura & D’Antonal 

IS 


(I2005I). Ivan d en Hoek & Groenewegenl d 19971 hereafter VHK), 
Marigd 1^0 iL hereafter MAR) , iGavilan, Buell. & Molla ll2005l); 
Gaydan. Molla. & Buelll l l200d . hereafter GAV), and Karakaj 
l 2010l hereafter KAR). A consequence of an ever-improving 
knowledge of AGB physics, is the reduction in the differences in 
the published yield s, in particular for the CN O elements. Thus, 
the work done by IStancliffe & Jeffe^ (|2007|) . centered on the 
mass loss rates, shows that changes in the yield by up to ~80% 
can result, but only for c ertain isotopes. On the other hand, 
IVentura & D’Antonal ( l2009l) focus their efforts on calculating new 
yields with significant l y impr oved values of the opacity. Finally, 
ICampbell & Lattanzid j200^ and iGil-Pons et al.l j2013h devote 
their work to the case of extremely metal poor stars, whose fi¬ 
nal evolutionary characteristics are not well known at the present. 
Apart from the AGB evolution, other authors have emphasised the 
importance of the nuclear reactions and their asso ciated numeri¬ 
cal pa rameters; this is the case for KAR an d also |Cristallp_et_aL 


(2009). The former re-derived the yields of iKarakas & Lattanzio 


1 20071) with new values of proton capture. The main differences 
with previous works reside in the yields of '^F, ^^Na and neutron 
rich isotopes^ Nevertheless, the CNO yields do not change signif¬ 
icantly amongst these works. 

Other imp ortant sets of yields available in the literature, such 
as lSiessl ( l20101) . have not been incorporated into our analysis. These 
authors provide yields for super-AGB stars with masses in the range 
7.5-10.5 Mq and metallicities between Z=le-4 and 0.04. The use 
of these tables implies a change in the mass at which one star 
explodes as a SN-II, nisN-n, and, more importantly, it introduces 
a third stellar mass range, instead of the two currently used (for 
low-l-intermediate and high-mass stars). We prefer for simplicity to 
adopt a constant, metallicity-independent, value of msN-u = 8Mq, 
rather than introduce an additional free parameter. We will explore 
the influence of this sort of metallicity-dependent SN-II mass limit, 
coupled with extant super-AGB yields, in the next future. 

Concerning the IMF, it is still matter of discussion if it is 
constant for all type of galaxies or if there are differences with 
environment, dependencies upon galactic stellar mass or metal- 
licity, or on the local star formation rate. Many recent works 


the regions (e.g. [Bekkil. l2013l: iGonrov et al.L l2013l: iDoncke et al 

2OI3I: Ferreras et all 20131: Geha et al,, 20131: Lasker et akL 

2013 

Me William. Wallerstein. & Mottinil 1201 3l: ISmith & Lucevl 

2013; 

Weidneretall 20131 and references therein), implying m 

most 


cases that the IMF might also be variable with tim e, but with dis¬ 
agreement among their results. ICalura et ^ l l2010t) used in a GEM 
an IMF which depends on the embedded cluster mass function, re¬ 
sulting in an IMF variable with time, as a function of the SFR. They 
conclude that the best fit to the solar neighbourhood data occurs 


* IChieffi & Limongl l2013l) also give new stellar yields but only for solar 
metallicity stars and for this reason they are not used here. 

^ The problems of and ^^ Na over-production f or AGB yields are out¬ 
lined in iRenda et al.l i2004l) and iFenner et al.l i2006l) . respectively. 


with an IMF resemb ling the standard one. lAndrews et al.l ( 1201 3h : 
IPeacock et al.l ( l2014h also support an invariant IMF for all types 
of systems. Regardless of these issues of invariance, the classical 
functional f orms for t he IMF employed in the lite r ature, includ¬ 
ing th o se of Salnetei (Il955h: Miller & Scald ( 197^ Ferrini et ^ 


(119901) : iKrounal j2002l) : IChabrieil ( l2003h and iMaschbergeJ ( 120131 
hereinafter SAL, MIL, FER, KRO, CHA, and MAS, respectively), 
whilst broadly similar, are quantitatively different from each other. 
In this work we will use these six forms, under the assumption that 
they are invariant with time. 

There are numerous CEMs in the literature, with important 
differences in their results, even for the case of the MWG for 
which the observational data sets are numerous. In these works, 
the selection of the best model, and the corresponding free pa¬ 
rameters, such as the star formation rate efficiency and/or infall 
rate, is performed, for any galaxy, comparing their observational 
data with a GEM built using a set of stellar yields with a given 

IMF (lGibsonLri997l:lGavilan. Buell, & Molli.l2005l : [Romano et~^ 

l20ld : Garigi & PeimbertL 201lh . Then, if observations cannot be 
well-reproduced, it can be claimed that an alternate set of y ields 
or IMF might be necessary dHemandez-Martmez et aI.Ll20I ih . Al¬ 
ternatively, it is possible to compare data with models computed 
using different IMFs to see which of these func t ions a re valid, 
without changing the stellar y i elds: iRomano et aP ( l2005l) did just 
that, concluding that iKroupal ( l200ll) . GHA. and MIL are better 
at reproducing the empirical data, than SAL, or IScald ( Il998h : 
IVincenzo, Matteucci, Belfiore & Maiolinoj ( l2015l) analyze the in¬ 
tegrated yields comparing results from different IMFs. However, 
the abundances within a galaxy or region therein, with a given star 
formation history, may be very different if another combination of 
IMF + stellar yields were to be used. Both ingredients are equally 
important to define the elemental abundances in a region and the 
corresponding temporal evolution. 

In this work, we make use o f the multiph a se chemical evo- 
lution model originally appli ed in iFerrini et alj d 19921 1 19941) and 
(Molla. Ferrini. & Diaa ( 119961) to the solar region, th e Galactic disc, 
and to other external spiral disks, respectively. In iMolla & DiazI 
(I2OO5L hereafter MD05), a large grid of models for a set of 440 
theoretical galaxies was generated. In that work the IMF was taken 
from FER and the stellar yields were from WOW and GAV. In ad¬ 
dition , the yields from Type la supernovae (SN-Ia) dlwamoto et al.L 
I1999I) w ere included along with the S N-Ia rate time distribu tion 
given bv lRuiz-Lapuente et alj dTOOOl) . In iGavichia et al.l d2014l) . we 
also used a similar model to that of MD05, applied to the MWG, 
modified to include bar-driven gas inflows, which has the effect of 
changing the SFR radial profile without significantly modifying the 
elemental abundance pattern. 

Our objective in this new work is to compute chemical evo¬ 
lution models for the MWG, with the same framework, total mass, 
molecular cloud and star formation efficiencies, and infall prescrip¬ 
tions for all of them, but with different combinations of stellar 
yields for massive stars (6 sets), low -1- intermediate mass stars (4 
sets), and IMFs (6 functions), thus resulting in a final grid of 144 
models. Our aim is to identify which is the best combination able 
to reproduce simultaneously the greatest number of observational 
constraints, mainly those pertaining to the radial distributions of 
gas, stars, and elemental abundances, and to the evolution of the 
solar region. 

The stellar yields and IMFs employed in our analysis are out¬ 
lined in Section 2. The chemical evolution model is presented in 
Section 3, along with the results of the 144 models. The selection 
of the best models is in Section 4, making use of a approach. 
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after comparison with the observational data (which are provided 
in Appendix ^At . Section 5 is devoted to our conclusions. 


2 INGREDIENTS: STELLAR YIELDS, INITIAL MASS 
EUNCTION, AND INTEGRATED YIELDS 


2.1 Stellar yield sets 


The stellar yield ^, (m) is defined as the fraction of the initial mass m 
of a s tar ejected in the form of freshly synthesised element i jPagel 
1 20091) 


qi(m) = 


ie,new,i 

m 


( 1 ) 


and is related to the total mass of this element i, niejejim), ejected 
by the star throughout its evolution (including pre-SN stellar winds) 
and death, via 


m^je,i(m) = m qi(m) + {m- mrc„,) A,-o, (2) 

where is the mass of the stellar remnant and A,,o is the abun¬ 
dance of the element i initially present in the star. 

Stellar yields are calculated by the stellar evolution commu¬ 
nity by coupling the evolution of the interior stellar structure with 
the relevant associated nuclear reactions. Such calculations provide 
the mass of each element produced and ejected to the ISM by stars 
of different masses throughout their lifetime. In chemical evolu¬ 
tion models, the stellar yields are usually divided into two ranges 
of stellar masses: 1) Low and intermediate mass stars, which in¬ 
clude those stars with masses m < 8 Mq; and 2) Massive stars, 
with m > 8 Mq, assuming that this is the minimum mass for stars 
which end their lives as SN-II. 


2.1.1 Low and intermediate mass stellar yields 

The main contribution from low and intermediate mass stars to the 
chemical enrichment is done during the Asymptotic Giant Branch 
(AGB) phase, where the mass-loss, thermal pulses. Third Dredge 
Up (TDU) events, and Hot Bottom Burning (HBB) are taking place. 
The first metallicity-depen dent yields used in GEM were those 
from iRenzini & VM ) ll98lh . In retrospect, the low mass loss rate 
adopted by the authors led to the need for a very large number of 
thermal pulses, to ensure reasonable remnant masses; the conse¬ 
quence of spending such a long time in the AGB phase was that 
almost of the ’^C was transformed into ''^N. 

As our knowledge of stellar evolution improved, newer yields 
were released with more accurate mass loss prescriptions, TDU 
events, and HBB. This is the case for the compilation of VHK, 
whose yields span a wide range of masses and metallicities (see 
Table[TJ, although still with very significant nitrogen production by 
stars with m > 4Mq. Later, armed with new stellar prescriptions, 
MAR calculated stellar yields for stars of masses between 1 and 
5 Mq. In her work, she suggested that stars with masses greater 
than 5 Mq end their lives as SN-II, thus only stars between 3 and 
5 Mq contributed to the nitrogen production. The final result was a 

sma ll excess in and a paucity of ''^nI^ _ 

iGavilan. Buell, & Moll^ ) l2005t) and iGavilan. Molla. & Buelll 


Table 1. Characteristics of the low and intermediate mass stellar yields used 
in this work. 


Set Name 

Z 

Mass Range 
(Mo) 

Yield 

Format 

Solar 

Abundances 

VHK 

0.001 

0.004 

0.008 

0.020 

0.8-8 

qi(m) 

AG89 

MAR 

0.004 

0.008 

0.020 

0.8-5 

mqi{m) 

GA91 

GAV 

0.0126 

0.0159 

0.0200 

0.0250 

0.0317 

0.8-8 

qiim) and nicjejim) 

GS98 

KAR 

0.0001 

0.004 

0.008 

0.020 

1-6 

qi{m) and rriejejiin) 

AG89 


AG89: [Anders & Grevess^ il989l): GA91: iGrevesse & AndersI jl99ll): 

GS98:[^ 


(120061 hereafter GAV) published new yields for low and intermedi¬ 
ate mass stars, with masses up to 8 Mq and a range of metallicities 
(see Table [T}. The main point of their work was the treatment of 
*^C and *'*N, concluding that a great amount of *^C in the ISM was 
ejected by intermediate stars, leading to '“^N yields not as great as 
VHK, nor as low as MAR, and reproducing well the observational 
constraints related to the time evolution of the elemental and rela¬ 
tive abundances of C, N and O, throughout the disc and halo0 

We use the stellar yields from VHK, MAR, GAV, and KAR 
for low and intermediate mass sta rs. Other excellent, mor e recent 
sets, such as ICristallo et alj ( 1201 ih or iLagarde et alj ( 1201 ih are less 
useful for our purpose here in that they either do not yet provide the 
full mass spectrum (the former compilation) or the CNO elements 
needed for our current work (the latter compilation). 

From these sets, VHK and MAR give their results as a frac¬ 
tion of stellar mass, g,(m), and as mass, mqi{m), respectively, while 
GAV and KAR give both, (net) stellar yields, mqi(m) ,and (total) 
ejected masses mgj^jim) (see Eq. 1). The relationship between both 
quantities depends on the initial abundances A,,o, usually assumed 
to be scaled to the solar ones for each value of the total abundance 
Z 


A,-,c 


Z 

■■ A,q —, 

'°Zg 


(3) 


while H and He abundances are take to be linear functions with Z: 


H 


He 


Up 


Zb ’ 


Hep + 


Hep - Hep ^ 
2q 


(4) 

(5) 


Thus, the initial abundances of each element, A,,o depend on 
the total abundance Z and also on the assumed solar values, A, q 
and f or the case of the H a nd He on the primordial values, Hp and 
Hep, djimenez et al.Ll2003l) as well. Solar abundances used are dif¬ 
ferent for each set, as specified in Table [T] the range of masses and 
metallicities are also listed there. We have interpolated the tables 


^ The impact of AGB yield selection, including iRenzini & Volil jl98lh , 
VHK, and MA R yields, as applied to CE M models of the Milky Way halo 
was explored bv iGibson & Mould jl997b . 


^ The evolution N/0 and its relation to 0/H is beyond the scope of this 
wo rk, but forms the ba sis of studies such as iGavilan. Molla. & Buelll j20Qq) 
an d lMollaetalJd2006l) . 
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Figure 1. Stellar yields for low and intermediate star. Each row shows an element. He, *^C, *^C, N, and O, from top to bottom, and each column refers to a 
different yield set as labelled. The number in each panel is the factor used to multiply the yields to plot all of them on a similar' y-axis scale. In each panel, the 
coloured lines represent different metallicities as labelled in the panel from VHK. 


given by these authors to obtain the ejected masses for the same 7 
metallicities: Z = 0.0, 0.0001, 0.0004, 0.004, 0.008, 0.02, and 0.05. 
We have also normalised the four sets, calculating comparable stel¬ 
lar yields for each. Tablej^gives these results for the suite of 
low and intermediate mass star yields employed here. 

The inferred solar abundances have (in large part) reduced 
from AG89 t o the most recent values, such as those from 
lAsplund et al.l ( l2009l) . Since the stellar evolution models employed 
here were constructed with the classic solar abundances, the yields 
must be used assuming that stars have those abundances. However, 
when analysi ng our results for the s olar region, we will use the most 
recent values ( lAsDlund et aUl200^ . 

To compare the different sets, we plot in Fig. [T] the stellar 
yields, as a function of the stellar mass, for He, *^C, *'^C, N 
and O. Although all sets show a broad similarity for each element, 
differences arises when we observe their behaviour in detail. As a 
generic result, MAR differs fhe most from the others, with a larger 
production for all elements and also a stronger dependence on Z, 
while VHK shows the smallest values. This is clear in the He pan¬ 
els, as for *^C, for which all sets show a maximum around 4 Mq 
and where MAR produces double the quantity of than KAR 


or GAV. For '^C, one can see a strong mass-dependence, with an 
abrupt increase for stellar masses only near 3 Mq. 

The stellar yield of N for these low and intermediate mass stars 
is very important since most of N proceeds from this stellar range 
and because a large contribution of the produced N is primary (NP): 
that is, independent of the original metallicity of the star. This NP is 
created in the HBB process, which needs a minimum core mass to 
initiate, as it occurs with the primary component of the *^C, as well. 
N appears for stellar masses around 4 Mq and when it appears, 
the consequently decreases. The behaviour is similar for all 
sets, except for MAR, which does not show, unlike the others, the 
increase at the highest mass. For O, the stellar yield is essentially 
negative (and very low in an absolute sense) for the entire mass 
range; only MAR shows positive values^ This negative yield will 
have consequences when the total integrated yield for oxygen is 
used, since the number of stars in this mass range is very high, 
compared with the number of the massive ones which produce the 
bulk of the oxygen. 


^ See also lPignatari et afl l2013l) for further recent evidence for modest pos¬ 
itive O stellar yields in this mass range. 
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Table 2. Stellar yields qiim) for our Urn sets. This is an example of the results for VHK and Z=0.02. The complete Tables 2a to 2d for VHK, MAR, GAV 
and KAR for the seven metallicities, are provided in the electronic edition. We give for each stellar yield set, the metallicity Z, and stellar mass, m, the stellai* 
yields, for elements as labelled, the remnant mass, nirem^ and, in the last two columns, the secondary contributions of and . 


z 

m 

Mq 

H 

D 

'He 

■’He 

‘'C 

"C 

N 

0 

fTlfgffi 

Mq 

"Cs 


0.02 

3.00 

-3.86E-02 

-3.81E-05 

-2.30E-05 

3.28E-02 

5.25E-03 

6.40E-05 

1.65E-03 

4.11E-04 

0.62 

4.15E-06 

1.37E-04 

0.02 

3.50 

-3.84E-02 

-3.83E-05 

-2.31E-05 

3.26E-02 

5.01E-03 

6.66E-05 

1.72E-03 

3.16E-04 

0.71 

3.47E-06 

1.15E-04 

0.02 

4.00 

-3.75E-02 

-3.85E-05 

-2.33E-05 

3.11E-02 

4.77E-03 

6.83E-05 

1.78E-03 

1.98E-04 

0.79 

2.79E-06 

9.47E-05 

0.02 

4.50 

-3.78E-02 

-3.89E-05 

-2.35E-05 

3.18E-02 

2.08E-03 

1.08E-04 

4.82E-03 

1.58E-04 

0.85 

2.82E-06 

9.89E-05 

0.02 

5.00 

-3.85E-02 

-3.92E-05 

-2.37E-05 

3.09E-02 

-5.94E-04 

1.40E-04 

7.24E-03 

1.42E-04 

0.92 

2.71E-06 

1.04E-04 



m (Mg) 

Figure 2. Ratio of primary to total stellar yield for and N: Top panels: 
Np/N; Bottom panels: *^Cp/*^C, for Urn sets by KAR, GAV, MAR, and 
VHK for dilferent metallicities, as labelled in Fig.[T] 


The contribution of the primary N in galaxies leads to the clas¬ 
sical relationship between N/O and O/H, in which a clear corre¬ 
lation for high metallicity/bright massive galaxies exists, but es¬ 
sentially none for low metallicity/low mass systems. Both con¬ 
tributions are separately given, or easily obtained, for GAV and 
MAR stellar yields, but for VHK and KAR it was necessary to 
calculate the primary contribu tion by the method described in 
iGavilan, Molli. & Buelll bOOhh . The ratio NP/N is shown in Fig.|^ 
for the four sets of low and intermediate stars used in this work. 
Obviously, the ratio NP/N is unity for Z = 0 and decreases as NS 
(secondary nitrogen) increases with Z. 


2.1.2 Massive star stellar yields 

The generation of massive star yields in the literature show the de¬ 
ployment of a range of evolutionary codes with different assump¬ 
tions regarding stellar micro-physics, including opacities and nu¬ 
clear reaction rates, and/or macro-physics, such as mixing or mass 
loss prescriptions. The NuGRID collaboration jPignatari et al.L 
l2013h has been established to rectify this heterogeneous situation, 
by employing an entirely homogeneous micro- and macro-physics 
approach across the full mass and metallicity spectrum (from low- 
to high-mass stars). However, at the time of pursuing this work, 
the only yields available publicly are for masses in the range [1.5- 
5]Mq and [15-60] Mq for Z=0.02 (for Z=0.01, the massive star 
range reduces to [15-25] Mq), without including the super-AGB 
phase and the SN-Ia stellar yields. Therefore, while the release of 


Table 3. Characteristics of the mas sets used in this work. 


Set Name 

Z 

Mass Range 
(Mo) 

Mass Loss 

Solar 

Abundances 

WOW 

0.000 

210-^ 

210-^ 

0.002 

0.02 

13^0 

N 

AG89 

PCB 

0.0004 

0.004 

0.008 

0.020 

0.050 

11-120 

Y 

AG89 

CLI 

0.0000 

io-« 

10-“* 

0.001 

0.020 

13-35 

N 

AG89 

KOB 

0.000 

0.001 

0.004 

0.020 

13^0 

Y 

AG89 

HEG 

0.000 

0.020 

10-100 

12-120 

Y 

LO03 

LIM 

0.000 

0.020 

13-80 

11-120 

Y 

AG89 


AG89:lAnders & Grevessj jl989l): L003:lLoddersll?003h 


this full grid is eagerly anticipated, it is premature to adopt it for 
these CEMs. 

Other well-known stellar yields are those that include a 
treatment of stellar rotation, the internal mixing and struc¬ 
tural changes resulting from which can lead to appreciable 
chan ges in the yields o f cert a in elem en ts, in particular nitro¬ 
gen I Kfevnet & Maedeii |2000L l2002al lbL IChiappini et al.L l2006l : 
iHirsctii . 2007h . A rich literature now exists which examines 
the role of this stella r rotation on stellar nucleosynthesis, al¬ 
though most of them jEkstr dm et al. L l2008l: iM evnet et Tl 201C 


I Yoon, Dierks. & Langer . 20121 : Chatzopoulos & Wheelei . 2012 
have emphasised the impact on very low-metallicity or Popula¬ 
tion III models. The lack of an available, fully self-consistent, grid 
of models spanning a range of mass and metallicity (up to solar) 
has somewhat restricted their application for chemical evolution 
studies. The precise treatment of rotational mixing, with veloci¬ 
ties varying from 60 to 800 km“‘ depending upon the authors and 
codes involved, remains a matter of debate. From a chemical evo¬ 
lution modeling perspective, the adoption of a given rotation veloc¬ 
ity (and its mass and metallicity dependence) implies an additional 
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free parameter, increasing the yield options dramatically^. Whilst 
acknowledging the importance of this issue, we feel it premature 
to proceed with a detailed comparison of the rotationally-mixed 
yields, until a fully self-consistent grid is available and calibrated 
unequivocally with empirical constraints. 

Thus, we use the extant compilations in the literature, in¬ 
cluding the sets from WOW, PCB, CLI, HEG, KOB, and LIM. 
In all cases, authors give their results as total ejected ma sses and 
most use solar abundances fromlAnders & Grevessa ( Il989h . except 
HEG who use iLodderd ( l2003h . This is given in column 5 of Ta¬ 
ble [3 where characteristics of the yields for each set are shown: 
the metallicities, the stellar mass range, and the inclusion or not of 
mass loss. 

As noted in Table [S] yields from WOW and CLI do not 
take into account pre-SN stellar winds and their corresponding 
mass loss, a fact particularly important for the most massive stars 
(m > 30 Mq), and give yields for an upper mass limited to 40 Mq. 
To extrapolate these yields up to 100 Mq is problematic, since mass 
loss is known to be substantial for these most massive stars. Ex¬ 
trapolation without the inclusion of mass loss would result in an 
integrated yield significantly higher than it should be. KOB, PCB, 
HEG, and LIM take into account a treatment of mass loss for these 
massive stars. In the three last sets, yields are provided up to 100- 
120 Mq, while KOB gives their results up 40 Mq. Therefore, we 
give in Table|4]the stellar yields up to 40 Mq for WOW, CLI, and 
KOB, and up to 100 Mq only for PCB, HEG, and LIM. Graphi¬ 
cally, we do the same in Figs and |4] in three panels (or columns 
of panels) at the left and at the right, respectively 

The stellar yields for these sets are given in Table |4] We com¬ 
pare, for the same elements as in Fig. [T] the different sets of yields 
of massive stars in Fig. [3] There, we show the results for each el¬ 
ement in a row, for each massive star yield set in a different col¬ 
umn, as labelled at the top of the figure. As a generic result, we see 
a very dilferent behaviour amongst the sets which are calculated 
without taking into account the existence of mass loss by stellar 
winds before the explosion of supernova (e.g., those of WOW and 
CLI), and those which do include these stellar winds (e.g., those of 
PCB, HEG, and LIM). Curiously, KOB is more similar to the first 
ones although formally this set forms part of the latter. The first sets 
are in the left-most columns, while the other four are in the right¬ 
most panels. The ones on the left show a weaker dependence on 
Z than the ones on the right, which is to be expected, as the mass 
loss is assumed to be dependent upon metallicity. In Fig. [3] we plot 
the yield of all elements multiplied by a factor as labelled in the 
WOW panel, in order to compare all of them on a similar scale. 
In all cases, H is negative with values in the range [-0.1,-0.5] de¬ 
pending on the stellar mass, on the metallicity, and on the authors. 
In LIM there is a strong variation around 30 Mq, and it becomes 
positive for the highest abundance. The behaviour for HEG shows 
quite abrupt changes with mass due to the Z = 0 set. 

Since He is produced directly from H, its behavior is comple¬ 
mentary to that of H, increasing when H decreases, although the 

® Besides that, it is not entirely sure to what degree, or if, this 
rotation is necessary for r eprodu cing observations. For example, 
Takahashi, Umeda, & Yoshidal ^2014^ have realized rotating and non¬ 
rotating models for Z = 0 for stellar masses between 12 and 140 Mq. Com¬ 
paring these models with the three most Fe-deficient stars in the Galaxy, 
they find that abundances for one of them are well-reproduced by 50-80 
Mq non-rotating models, the second one is equally well-fitted with non ro¬ 
tating or rotating 15-40 Mq models, and only one of them might require 
rotating 30-40 Mq models 



10 20 30 40 20 40 60 80 100 

m (Mq) m (Mg) 

Figure 4. Stellar yield for n-elements a Ne Mg Si S Ca 
for mas sets for different metallicities. Lines are coded as in the previous 
figures. 

absolute value is smaller than that, since a certain quantity is nec¬ 
essary to create the other elements. The ejected mass of He and 
is higher in the case of PCB, HEG, and LIM than those of WOW, 
CLI, and KOB. When we compare the same mass range we see 
that for m < 40 Mq more He, C, and N is ejected, while O is pro¬ 
duced in a smaller quantity when the mass loss by stellar winds is 
included. 

For the elements beyond O, the yields are shown in Fig. |4] 
for the six different sets of massive stellar yields. Here, we add 
the yields for elements Ne, Mg, Si, S, and Ca and represent this 
Qf-yield as a function of the stellar mass for each yield set. As for 
O, a different behaviour arises between the yields calculated tak¬ 
ing into account the stellar winds (right panels) and those which do 
not (left panels); the former show a maximum around 20-30 Mq. 
KOB shows a behaviour between both, similar to those without 
mass loss, but also indicating a slight maximum near 40 Mq. 

2.2 The Initial mass function 

We are building upon the Galactic model outlined in MD05, but in¬ 
stead of simply using the PER IMF (as we did in that work), we now 
employ a range of functional forms for the IMF, as well as the var¬ 
ious stellar yield data sets described in §2.1. The IMFs adopted are 
from SAL, MIL, FER, KRO, CHA, and MAS, as shown in Fig. [5] 
where differences amongst them appear readily. We assumed the 
IMF to be invariant with time and metallicity. 

The functional forms for the adopted IMFs are: 


f){m)sAL 

,.^-2.35 
= m , 


(6) 

4'dn)MIL 

(logm-Fl.02)2 

= e 2 0.68- , 


(V) 

(p(rn)fER 

_ JQ-V'^ ’^^+logn 

(1.92+2.07 log m)/„^1.52 

(8) 




0.15 < ot/Mq < 0.08 


^(tn)KRo 

= . 

0.08m-‘-3 

0.04m-2-3 

0.08 < ot/Mq < 0.50 

0.50 < m/MQ < 1 

(9) 



0.04m-2'’ 

m/MQ > 1 
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Figure 3. Stellar yield for H, He, *^C, *^C, N and O for massive stars from sets by WOW, CLI, KOB, PCB, HEG, and LIM, for different metallicities coded 
with colours as in Fig.[T] 


Table 4. Stellar yields qi(m) from different massive stars mas sets. The values for WOW and Z=0.02 are given here as an example. The complete set of tables 
for WOW, PCB, CLI, KOB, HEG, and LIM, for seven metallicities (Z=0, 0.0001, 0.0004, 0.004, 0.008, 0.02, and 0.05) are provided in the electronic edition. 


z 

m 

Mq 

H 

D 

"He 

“He 

^C 

mo 

0.02 

22.00 

-0.240E+00 

-0.436E-04 

-0.701E-06 

0.908E-01 

0.823E-02 

O.lOOE+00 

0.02 

25.00 

-0.270E+00 

-0.440E-04 

-0.144E-05 

0.925E-01 

O.lOlE-01 

0.122E+00 

0.02 

30.00 

-0.309E+00 

-0.449E-04 

-0.319E-05 

0.886E-01 

0.686E-02 

0.155E+00 

0.02 

35.00 

-0.335E+00 

-0.452E-04 

-0.472E-05 

0.801E-01 

0.632E-02 

0.158E+00 

0.02 

40.00 

-0.331E+00 

-0.415E-04 

-0.570E-05 

0.867E-01 

0.648E-02 

0.143E+00 


""Ne 

"“Mg 

^^Si 

""S 

^'’Ca 

"'’Fe 

nirem 

"Cs 


0.319E-02 

0.356E-03 

0.157E-01 

0.723E-02 

0.423E-03 

0.841E-03 

2.02 

0.375E-04 

0.206E-02 

0.158E-01 

-0.113E-02 

0.121E-01 

0.504E-02 

0.301E-03 

0.821E-03 

2.07 

0.320E-04 

0.216E-02 

0.144E-01 

0.782E-02 

O.lOOE-01 

0.272E-02 

0.692E-04 

0.787E-03 

1.94 

0.265E-04 

0.244E-02 

0.254E-01 

0.661E-02 

0.271E-02 

-0.132E-03 

-0.335E-03 

0.738E-03 

2.03 

0.205E-04 

0.254E-02 

0.310E-01 

0.429E-02 

0.916E-03 

-0.274E-03 

-0.308E-03 

0.658E-03 

5.45 

0.161E-04 

0.257E-02 
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log m (Mg) 


ntfliip 

N, = I A^{m)dm. (16) 

Jmi„ 

Our initial plan was to use the same lower =0.15 Mq) 
and upper (m,,,, = 100 Mq) mass limits for each CEM; how¬ 
ever, as noted previously, some yield compilations are restricted 
to <40 Mq. As such, we have computed the number of stars for 
each IMF for these two values of m,,,, (see Table and in the 
next section, models have been computed for each combination of 
IMF+massive stars with a different m„p following the set of mas¬ 
sive stars used. 


3 CHEMICAL EVOLUTION MODELS 


Figure 5. The initial mass functions used in this work as log m^(m) by SAL, 
MIL, FER, KRO, CHA, and MAS. 


Table 5. Number of stars for the adopted IMFs for a stellar mass of 10 "^Mq . 
We show the normalisation constant A, the total number of stars, N,, the 
number of stars with mass smaller than IMq, Niow, the number of low 
and intermediate mass stars, with 4Mq < m < 8Mq, Num, the number of 
massive stars with m > 8Mq, which will he SN-II, Nsn, and the number 
of stars more massive than 20 Mq, 


IMF 

A 

N, 

^iow 

Nlim 

Nsn 

^mas 

SAL 

2090 

m,p 

20200 

= 40Mo 
19970 

145 

83 

17 

MIL 

191 

13387 

13093 

214 

79 

9 

FER 

22000 

18924 

18792 

98 

34 

5 

KRO 

80830 

17100 

16925 

124 

51 

8 

CHA 

148210 

11125 

10780 

215 

129 

26 

MAS 

13110 

13133 

12800 

206 

125 

25 

SAL 

2000 

ntup ■■ 

19715 

= 100 Mq 
19490 

138 

86 

23 

MIL 

189 

13329 

13037 

212 

80 

11 

FER 

21869 

18837 

18704 

98 

35 

6 

KRO 

79458 

16932 

16756 

123 

54 

10 

CHA 

137808 

10386 

10055 

200 

131 

36 

MAS 

12222 

12319 

11999 

192 

128 

35 


^{m)MAS 


0.043m 


AA 


f^char 


(k)gHi-bQ.657)^ 


1 + 


^char 


0.15 < hiIMq < 1 

(10) 

1 < mIMQ < 100, 




• 

(11) 


(12) 


where: rrichar = 0.2Mq, 

, I m 

G(m) = I 1 + 

AA = 


^char 

(1-a) (1-/3) 


and 

1 


J r>m 


m,har G{m„p) - G(mia„) 

As usual, the total mass in stars is normalised to 1 Mq 
A m (j){m) dm = 1 Mq. 


(13) 

(14) 

(15) 


and in this way, the total number of stars, A, in a generation is 


3.1 Summary Description 

The chemical evo lution c ode used here is tha t described in MD05 
and IMoII j ( l2014h . and in iMolla et al.l ( 1201 Sl hereinafter MCGD), 
the latter in which we present a new updated grid of chemical evo¬ 
lution models for spiral, irregular, and low mass galaxies with some 
modifications in the input parameters over the ones from MD05. 

We assume a radial distribution of primordial gas in a 
spherical proto-halo falls onto the plane defining the disc0 The 
mass radial distribu tions are calculated from the prescriptions in 
ISalucci et alj ( l2007h . who give expressions to compute the halo 
density, virial radius, rotation curve, and final mass of the disc as 
functions of the virial mass, Myi^. We have calculated an initial mass 
distribution with a dynamical mass of 10*^ Mq and a maximum 
rotation velocity of Vrot = 177kms“'. The infall rate or collapse 
timescale in each radial region is chosen in such a way that the disk 
ends with a radial profile similar to that observed, by following the 
prescriptions from ishankar et all (l2006l) for the ratio Mdyn/Mdisk, 
at the end of the evolution for a time of 13.2 Gyr. This method gives 
as a result, for the chosen virial mass, a radial distribution of the fi¬ 
nal mass of the disk Mo{R) and also the collapse timescale radial 
distribution, t{R), necessary to obtain it. 

Our formalism for the SFR adopts two stages, first forming 
molecular clouds from the diffuse gas according to a Schmidt law 
with n = 1.5, and then second, forming stars from cloud-cloud col¬ 
lisions. Once choosing the total mass radial distribution, it is nec¬ 
essary to determine which are the best efficiencies to form molec¬ 
ular clouds and stars for this MWG-like galaxy. This has been per¬ 
formed in MCGD, comparing the time evolution of the region lo¬ 
cated at R=8 kpc and the radial distributions of gas, stars, and SFR 
with the present-time data. This comparison allowed us to select 
and fix the best efficiencies to reproduce the MWG disk data, which 
are the one used in this work. 

For this basic model we have computed all possible combina¬ 
tions of the six IMFs with the six mas sets and the four Jim yields 
described in ^ resulting in a total of 144 models for the MWG. 
In order to identify the best combination capable of reproducing 
the extant observations, we will now compare the results of these 
models with the observational data given in Appendix jAl where, 
furthermore, the empirical data has been binned in order map these 
onto our model bins. 


^ The code is inherently one-dimensional (in r), involving a thin disc and 
azimuthal symmetiy. 
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Figure 6. The radial distributions of the MWG disk surface densities for the 
sets of 144 models compared with the observational data as red dots with 
error bars: a) the surface density of diffuse gas, Lhi, b) molecular gas, 
c) stars, Z*, in units of Mq pc“^ for all of them, and d) SFR, "Lsfr^ in units of 
Mq pc“^ Gyr“^ All panels are given in logarithmic scale. Each colour-type 
of line indicates a different IMF with the same coding as in Fig. [5] 


Figure 7. The evolution of the Solar Neighbourhood for the set of 144 
models compared with observational data as red dots with eiTor bars, as 
obtained in § IaT] The large yellow dot represents the solar values, a) SFR 
(in MQyr^9 in logarithmic scale; b) [Fe/H]; c) [cr/Fe] vs [Fe/H]; d) The 
MDF. The coding of the lines is as in Fig. [6] 


3.2 Results for the Solar Vicinity and the MWG disc 

In Fig. 1^ we present, for the 144 models, the results concerning 
the state of the disc or the radial distributions at the present time 
(t = 13.2 Gyr) for gas, stars, and the SFR, compared with the data 
shown in Table lA^ For the SFR, panel d), we have artificially in¬ 
creased all values (models and data) by 1 dex, in order to plot them 
in the same scale as used in panel c). The results show a small 
dispersion around the data or around the mean values. These good 
results are the consequence of the infall rate and the star forma¬ 
tion rate efficiencies selected for the model to reproduce the MWG. 
In all cases the models’ dispersions are comparable to, or smaller 
than, the data uncertainties. These radial distributions are, as ex¬ 
pected, only slightly dependent upon the IMF, due to the different 
rate of ejected/retumed gas by (mostly massive) stars when they 
die. 

The evolution of the SFR, metallicity, [tt/Fe] as a function of 
[Fe/H], and the Metallicity Distribution Function (MDF) for the so¬ 
lar region of our 144 models are compared with the observational 
data given in Tables [ATI and IA2l in Fig.|3 It is quite evident that, 
even with the same input parameters and total mass for MWG, the 
resulting evolution is different for each model. In the case of the 
SFR, this is due to the IMF used, since the value of the mass locked 
in stars (1-R) changes with IMF. Thus, the evolution for FER shows 
the lowest SFR histories while MAS models show the highest, since 
the returned gas fraction of each stellar generation is the lowest 
for the FER models. Within each IME, each combination mas+Iim 
also shapes the results somewhat, but in this panel the yields have 
a smaller role than the IMF. In panel b), the results are the con¬ 
sequence of the SFR history, and, therefore, again FER is in the 
lowest part of the model locus while MAS is in the highest. Since 
Fe is produced mainly by the SN-Ia, the results depend more on 
IMF than on the stellar yields of massive stars; they do not depend 
upon the low- and intermediate-mass stars. 

In panel c) we show the classical plot of [a/Fe]-[Fe/H] for the 
solar region, where [Fe/H] is often taken as a proxy for time. This 
figure gives the differences in the ejection to the ISM of ff-elements, 
coming from massive stars, and from the Fe ejected mainly by SN- 


Ia, and also partially due to mas yields. Therefore, both IMF and 
massive stellar yields are playing a role here. It is evident that a ‘by 
eye’ inspection of these panels would suggest that the KRO, CHA, 
MIL, and SAL in our models reproduce better the data. When we 
use MAS, results fall above the data for all combinations of stellar 
yields, while our models using FER tend to lie below the obser¬ 
vations. This plot also gives an indication concerning the massive 
star yields + IME combinations which may be rejected: WOW is 
only valid when used with FER. In fact, WOW have already noted 
that their Ee yield is high and recommend it be divided by a factor 
of two in order to best reproduce the data with a GEM. This high 
Fe yield is compensated for when using FER, since the number 
of massive stars is small in this IMF compared with the others. In 
panel d) it is again evident that the IMF has an important effect on 
the MDF, with most of FER models at the left and MAS models at 
the right of the observations. Again our KRO, MIL, and some CHA 
models seem to fit better the observed MDF. 

Finally, we present the resulting radial distributions of C, N, 
and O for the whole set of 144 models in Fig. [8] compared with 
the binned data obtained in 3 A3 1 The slope of the radial abundance 
gradients does not depend, as expected, on the combination IMF+ 
stellar yield. The radial gradient is determined by the ratio between 
SFR and infall rate, 'V(t)lf, and it is basically independent of the 
IMF or the stellar yields. This is the reason why the radial gradient 
is basically the same for the 144 models, since we use the same SFR 
and infall history for all of them; that said, the absolute abundances 
change significantly, since, even using the same basic model, the 
combination stellar yields+IMF may change the absolute values of 
abundances in the disc by a factor of 100 for C, more than a factor 
of 10 for N, and a factor of 30 for O. Thus, the 144 models results 
show a dispersion clearly greater than the data and the comparison 
with data allows us to select the appropriate combination of yields 
and IMF. 

This is the most important result of this section: that it should 
be possible to select, on the basis of our GEM, which of these com¬ 
binations may be valid in reproducing the empirical data and which 
of them should be rejected. This is an important point to note as, 
in order to reproduce a given observation which appears to show a 
higher or lower metallicity than predicted by a model, a common 
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Table 6. Values of obtained from the fitting of models to each one of the data sets defined in Appendix [aI The entire table is presented in the electronic 
version. A portion is shown here for guidance regarding its form and content. 


Urn 

mas 

IMF 


- 1 - 

Ffr-e///! 

\alFe\ 

-7- 

^MDF 

^H1 

Fh, 

-7 — 


^CIH 

^NIH 

^OIH 

GAV 

CLI 

SAL 

11.768 

2.387 

70.848 

92.836 

1.623 

5.157 

2.517 

2.805 

51.846 

7.614 

9.563 

GAV 

CLI 

MIL 

17.216 

3.894 

1.338 

112.469 

1.803 

5.462 

1.227 

3.669 

69.440 

13.204 

3.302 

GAV 

CLI 

FER 

19.522 

34.283 

43.317 

158.274 

1.982 

5.070 

4.553 

2.593 

23.120 

29.962 

21.974 

GAV 

CLI 

KRO 

14.206 

6.894 

1.665 

50.123 

1.666 

5.068 

3.114 

2.743 

31.778 

5.876 

3.925 

GAV 

CLI 

CHA 

15.985 

6.429 

48.694 

167.581 

1.728 

5.411 

1.283 

3.520 

77.556 

21.434 

13.190 

GAV 

CLI 

MAS 

22.255 

18.741 

154.480 

364.015 

2.014 

5.719 

0.690 

4.135 

119.605 

58.290 

42.565 


Table 7. Confidence levels for the 8 best models. For each combination Iim+ mas + IMF, defined in columns 1, 2, and 3, the confidence levels obtained 
when fitting separately each data set of observations to our models, for columns 4 to 14. Column 15 is the combined likelihood, Pn, calculated using all 
observational sets. Column 16 is Pj, eliminating the disc properties (stars, gas, and SFR radial distributions). 


Urn 

mas 

IMF 

P't, 

P[Fe/H] 

P[a/Fe] 

Pmdf 

Phi 

F«2 

p. 

P'is 

Pc/H 

Pnih 

Pq/h 

Fit 

Pi 

VHK 

CLI 

KRO 

0.997 

0.913 

1.000 

0.502 

1.000 

1.000 

0.994 

1.000 

0.997 

0.999 

1.000 

0.930 

0.894 

MAR 

CLI 

KRO 

0.997 

0.972 

0.995 

0.770 

1.000 

1.000 

0.994 

1.000 

0.799 

0.443 

1.000 

0.885 

0.826 

GAV 

CLI 

KRO 

0.996 

0.881 

0.994 

0.267 

1.000 

1.000 

0.992 

1.000 

0.488 

1.000 

1.000 

0.819 

0.732 

KAR 

CLI 

KRO 

0.998 

0.993 

1.000 

0.108 

1.000 

1.000 

0.994 

1.000 

0.998 

1.000 

1.000 

0.816 

0.727 

MAR 

HEG 

MIL 

0.729 

0.968 

0.129 

1.000 

1.000 

1.000 

1.000 

1.000 

0.383 

1.000 

1.000 

0.737 

0.619 

MAR 

HEG 

KRO 

0.998 

0.729 

0.410 

0.379 

1.000 

1.000 

0.996 

1.000 

0.372 

0.523 

1.000 

0.707 

0.580 

VHK 

PCB 

FER 

0.950 

0.697 

0.970 

0.035 

1.000 

1.000 

0.983 

1.000 

0.533 

1.000 

1.000 

0.668 

0.532 

MAR 

KOB 

KRO 

0.996 

0.530 

0.825 

0.069 

1.000 

1.000 

0.993 

1.000 

0.995 

0.265 

0.999 

0.644 

0.501 
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Figure 8. The radial distributions of elemental abundances, as 12 + 
log(X/Fl): a) C, b) N and c) O, in the MWG disx for the sets of 144 models, 
compared with observational data as red dots with error bars. Coding of the 
lines is as for Fig.|3 


fall-back option is to invoke some mechanism(s) of mixing, enrich¬ 
ment, or dilution of abundances, to reconcile the discrepancy. As 
we show here, the alternate suggestion that the correct selection 
of IMF or/and stellar yields may be an easier way to achieve the 
desired abundance patterns should not be dismissed. 


4 THE SELECTION OF THE BEST MODELS 
4.1 The application of a technique 

The objective of this section is to find the best combination of IMF 
+ stellar yields able to reproduce the MWG data amongst the 144 
models computed and described in the above section. In order to 
do this, we use a classical x^ technique comparing the model re¬ 
sults and the corresponding observational data, such as those used 
in Fig. |6] |7] and [8] In Table |6] we give our results; for each 
model calculated with a combination of Urn set +mas set and IMF, 
we show the x^ obtained from the comparison of our models with 
the data for all observational sets we use. 

As said before and shown in Fig. [6] all models are equally 
good at fitting the radial distributions of both phases of gas, stars, 
and the SFR. We confirm this fact with the values of x^ for these 
quantitiel^ in Table Basically for all models they fall below the 
limits corresponding to 80% of confidence level; that is, models 
fulfill widely these constraints. Therefore, we analyse the fit of our 
models for the other 7 empirical datasets. 

We have assumed that each model is represented by ax^ dis¬ 
tribution, and calculated the corresponding likelihood, F,, or confi¬ 
dence level, (complement of the significance level a associatetj^ to 
each^^). The number of free parameters, NF = 3 in all cases, and 
the number of points for the fitting, variable for each data 

set i, give the number of degrees of freedom fc, = NF - The 
likelihood is calculated as 


* In all cases, the value at R = Okpc has not bern used in these calcula¬ 
tions, since the differences between data and models are large in this region, 
and thus our x^ values would be biased toward models with high densities 
in the inner disk, regardless of the quality of the agreement elsewhere in the 
disc. 

^ cr is the statistical significance, corresponding to a given giving the 
probability of rejecting the null hypothesis, given that it is true, the null 
hypothesis being that both sets (observations and model results) would rep¬ 
resent the same sample. 
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[Fe/H] 


Figure 9. Solai* vicinity evolution of the best 4 models shown in Table [7] 
a) 'PCO, b) [Fe/H](t), c) [o'/Fe] vs [Fe/H] and d) MDF Red and yellow dots 
have the same meaning as in Fig. [7] The orange dot-dashed line represents 
the model where Pmz>f=1 but P[Q./Fe] = 0- 


Pi = i- a(Xk. <x) = 




After computing these likelihood values, we see that the SFR 
and enrichment histories, much like the C, N and O abundances, 
may be easily reproduced with some combinations of IMF+ yields, 
showing low values forx^, and high likelihood Pi values. However, 
the relation [Q'/Fe]-[Fe/H] and the MDF are more difficult to fit, and 
thus constrain the selection of models able to reproduce simultane¬ 
ously all data sets. 

In order to choose the best models, we have computed the 
combined likelihood, Ps 


Ps 


s 

n^- 


V;=i,M 


(18) 


obtained as the geometrical average of the individual Pi previously 
calculated for each data set, and S is the number of used data sets. 
In this expression, we may assume that a good model is the one that 
simultaneously reproduces all data sets, including the ones pertain¬ 
ing to the present state of the disc, 'Lhi, 2., and 'Lsfr\ in that 
case, the number of datasets used is 5 = 11. Conversely, we could 
only use the 7 datasets shown in Fig.|7]and[8] that is, the observed 
SFR and enrichment histories, the relation [a/Fe]-[Fe/H], the MDF, 
and the radial profiles of C/H, N/H, and 0/H. Therefore, by max¬ 
imising the combined likelihood P^ or P-j, we are able to select the 
best models of our grid. We have computed both values Fn and Fy, 
and then, we have ordered our models by using the combined like¬ 
lihood Fy and taken the first 8 (which represents 5% of the total 
number of the calculated models), which have values Fy>~50% . 
The order of models using Fn is exactly the same for these mod¬ 
els, showing values Fn > 65%. We show these models in Table|7] 
Only four from our models present values higher than ~70% for 
the fit in the seven selected data sets (or > 80% in the entire set of 
observations) and all of them use CLI-l-KRO combinations. All the 
other models have Fy < 50% . 

These results allow us to constrain the models, reducing the 
valid ones to only 4-8 models, depending on the goodness we re¬ 
quire. However, we must take into account that this conclusion is 
mainly due to the MDF, which show values of very high com¬ 



Figure 10. Radial distributions of elemental abundances for the best 4 mod¬ 
els compared with data. The meaning of colours, symbols, and types of lines 
is the same as in Fig.|9] 


pared with most of the models. In fact, besides MAR+HEG+MIL, 
the 5"' in the table, there there is only one o ther model, corre¬ 
sponding to the combination MAR-l-HEG-l-SAL, which has a value 
PMOf=l- However, these two models have a P[a/Fc]=0.129 and 0, 
respectively, which implies they do not reproduce this relation at 
all. Actually there is only one model, the first one of the table, with 
F, > 50% for all columns. If we eliminate the MDF as a constraint 
for our models and calculate the equivalent P^ and Fjo, we find 11 
models satisfying this condition for the ten other columns. Eight 
of them, shown in Table [D using MIL or EER as IMF, are able 
to reproduce the six (or ten) remaining data sets within a confi¬ 
dence level Fj higher than 80% (or ~ 90% for F^). In fact, mod¬ 
els 1 to 4 in Table |7] showing Fn > 80%, increase to values 
Fio > 92%, when we don’t take into account the columns corre¬ 
sponding to the MDF. Therefore, models of this second table would 
also be valid, considering that many of the litera ture MDFs of the 
past d ecade should different maxima positions: ICasagrande et al.l 
(I2OIII1 found this m aximum at [FelH] ~ -0.0 5 , similarly 
to IChang. Hou. & FtJ ||20 0(]|): iLuc k & HeiteJ ( l2006h : iFuhrmann 
( I2OO8I) , while Kordopatis et al. I2O I5I) fi nd it near -0.2 dex, more 
in agr e ement wit h lAllende Prieto al.| ([2004|): Nordstrom et al.l 
( 1200^ : [Holmberg. Nordstrom. & AndersenI ( |2007|) . Although we 
have used some of these datasets to obtain a bin-averaged MDF, 
it is likely that the error bars associated with these data are higher 
than the pure statistical ones included in our^^ calculation. 

Summarising, our best models are combinations of CLI-KRO 
with any lim set. It is necessary to note that, given the possible un¬ 
certainties in the MDF, perhaps other combinations of stellar yields 
and IMF, as shown in Table [S] might succeed in reproducing the 
MWG data, mainly if other hypotheses pertaining to the evolution¬ 
ary scenario (infall rate or SFR) are assumed. 


4.2 Results for the best models 

Having selected our best models, we plot their results in the sub¬ 
sequent figures to compare with the observational data. In Fig. 
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Table 8. Confidence levels for eight other best models selected without using the MDF. 


Urn 

mas 

IMF 

P't, 

P[Fe/H] 

P[a/Fe] 

Pmdf 

Phi 

Phi 

p. 

P^. 

Pc/H 

Pn/h 

Pq/h 

Pio 

Pe 

KAR 

KOB 

MIL 

0.915 

0.998 

0.998 

0.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

0.991 

0.985 

KAR 

CLI 

MIL 

0.873 

0.939 

0.919 

0.000 

1.000 

1.000 

1.000 

1.000 

0.997 

1.000 

0.999 

0.972 

0.953 

MAR 

KOB 

MIL 

0.908 

0.993 

1.000 

0.000 

1.000 

1.000 

1.000 

1.000 

0.747 

1.000 

1.000 

0.961 

0.936 

KAR 

PCB 

FER 

0.958 

0.800 

0.977 

0.009 

1.000 

1.00 

0.982 

1.000 

0.880 

0.996 

1.000 

0.957 

0.932 

MAR 

CLI 

MIL 

0.856 

0.997 

0.999 

0.000 

1.000 

1.000 

1.000 

1.000 

0.465 

1.000 

1.000 

0.912 

0.857 

KAR 

KOB 

KRO 

0.997 

0.661 

0.976 

0.002 

1.000 

1.000 

0.992 

1.000 

0.532 

0.999 

1.000 

0.898 

0.836 

VHK 

PCB 

FER 

0.950 

0.697 

0.970 

0.035 

1.000 

1.000 

0.983 

1.000 

0.533 

1.000 

1.000 

0.897 

0.836 

VHK 

KOB 

KRO 

0.994 

0.310 

0.993 

0.016 

1.000 

1.000 

0.993 

1.000 

1.000 

0.999 

1.000 

0.887 

0.821 


we show the evolution with time of SFR, 'P(t), [Fe/H](t), [a/Fe]- 
[Fe/H], and the MDF. We have also drawn as an orange dot-dashed 
line the model MAR+HEG-fSAL which does not reproduce the re¬ 
lation [cr/Fe]-[Fe/FI] as said in the previous section. 

Finally, we show in Fig. [TO] the elemental abundances of a) 
C, b), N, and c) O with the same line coding that in the previous 
Fig.ia We see in the panel b) that model using VFIK shows the 
highest N abundances of the four models, just within the limit of the 
uncertainties, while using MAR, with G AV lying between the two 
and cl osest to the date, as also found in iGavilan, Molla, & Buelll 
( I2OO6I) . T he radial gradient o f O abundances obtain ed with recent 
data from iFtenrv et al.l ( 1201 Oh and lLuck et al.l ( 1201 ih gives an aver¬ 
aged value of -0.040 dexkpc^*; C data gives a radial gradient of 
-0.048 dexkpc^', similar to the one for O. For N we obtain a ra¬ 
dial gradient of -0.062 dexkpc^*, slightly steeper than the one for 
C and O. The four models show radial distributions which seem in 
fair agreement with the observed radial gradients. 


5 CONCLUSIONS 

• By using our multiphase chemical evolution code, we have 
calculated 144 models applied to the MWG, with the same basic 
hypotheses, but different combinations of four low and intermediate 
mass stellar yield sets, with six massive stellar yield sets, and six 
IMFs. 

• We have analysed the observational data corresponding to the 
temporal evolution for SFR and iron abundance, the relative abun¬ 
dance [or/Fe] as a function of [Fe/H], and the MDF for the solar 
region; further, we provided radial distributions of masses and ele¬ 
mental abundances at the present time for the Galactic disk, obtain¬ 
ing binned data sets averaged with different authors’ samples. 

• Using a classical technique, we compared the results of our 
144 models with the binned data points from the observational data. 

• Assuming that a good model is the one that simultaneously 
reproduces the observed SFR history, the [a/Fe]-[Fe/H] relation, 
the MDF, and the radial profiles of C/H, N/H, and O/H, we defined a 
geometrical averaged likelihood from the product of the individual 
confidence levels for these 7 quantities. 

• We find that the best 4 of our 144 models are able to reproduce 
all observational data sets with confidence levels P^ higher than 
~70%, and use combinations CLI-l-KRO with any Jim yields. It is 
necessary to take into account that, given the possible uncertainties 
in the MDF, maybe other different combinations of stellar yields 
and IMF might be equivalently good to reproduce the MWG data, 
mainly if other assumptions regarding the infall rate or SFR are 
used. 
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Table Al. Binned SFR and metallicity evolution for the solar vicinity 


Time 

(Gyr) 

log T Ai/ogT) 

(Moyr-‘) 

[Fe/H] 

A([Fe/H]) 

0 

-0.3143 

0.1667 

-0.788 

0.25 

1 

-0.0631 

0.0935 

-0.208 

0.26 

2 

-0.2668 

0.2510 

-0.225 

0.25 

3 

0.0165 

0.0507 

-0.184 

0.24 

4 

-0.0118 

0.1055 

-0.097 

0.25 

5 

0.0127 

0.0577 

-0.073 

0.25 

6 

-0.0651 

0.1568 

-0.040 

0.23 

7 

-0.1400 

0.1372 

0.017 

0.22 

8 

-0.0530 

0.1534 

-0.007 

0.22 

9 

-0.2535 

0.1797 

0.005 

0.22 

10 

-0.2420 

0.2441 

0.056 

0.21 

11 

-0.4513 

0.0809 

0.074 

0.21 

12 

-0.3280 

0.1581 

0.077 

0.21 

13 

-0.6039 

0.3078 

0.160 

0.23 

13.2 

-0.6676 

0.3759 

0.200 

0.23 


Table A2. Binned [Q'/Fe]-[Fe/H] relation and MDF for the solar vicinity 


[Fe/H] 

[ff/Fe] 

A[a/Fe] 

MDF 

Amdf 

-1.50 



0.039 

0.030 

-1.40 



0.062 

0.040 

-1.30 



0.063 

0.042 

-1.20 



0.013 

0.050 

-1.10 



0.072 

0.038 

-1.00 



0.129 

0.050 

-0.90 



0.127 

0.042 

-0.80 



0.190 

0.100 

-0.70 



0.147 

0.030 

-0.60 

0.268 

0.019 

0.411 

0.088 

-0.50 

0.205 

0.005 

0.502 

0.013 

-0.40 

0.176 

0.002 

0.886 

0.123 

-0.30 

0.134 

0.005 

1.440 

0.169 

-0.20 

0.071 

0.004 

1.900 

0.065 

-0.10 

0.025 

0.003 

1.870 

0.077 

0.00 

-0.021 

0.004 

1.500 

0.141 

0.10 

-0.051 

0.009 

0.971 

0.172 

0.20 

-0.097 

0.008 

0.769 

0.154 

0.30 

-0.115 

0.009 

0.239 

0.057 

0.40 

-0.140 

0.011 

0.076 

0.017 

0.50 

-0.122 

0.021 



0.60 

-0.145 

0.012 



0.70 

-0.200 

0.000 




APPENDIX A: OBSERVATIONAL DATA 

The observational data against which our CEMs are compared are 
now outlined. These include the solar neighbourhood’s temporal 
evolution, in addition to the present state of the disc, including the 
radial distributions of surface densities for stars, gas, and star for¬ 
mation rate, and elemental abundances of C, N, and O. Other data, 
such as [X/Ee], are usually represented as a function of [0/H] or 
[Ee/H], with the latter typically being employed as a proxy for time. 
Thus, we have also used the [ff/Re] - [Re/H] data of the solar vicin¬ 
ity to compare with our models. 

Al The Solar Vicinity 

Ror the solar vicinities of our CEMs, we compare with extant ob¬ 
servations pertaining to the time evolution of the SRR and the 
age-m etalli city relation. The SER ev olution is taken from iTwarod 
( 1 19801) and [Rocha-Pinto et alj l l200(t) . In both cases, the data show 
a maximum around 8-10 Gyr ago, that is, the onset of the SER 
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Figure A l. Solar neighb ourh ood data: a) Star formatio n history T(/) with 
data from lTwarod il98(]h and iRocha-Pinto et alJ i2Q0Qh as blue and green 


from iTwaro 

3 fi 98 cl); 

Edvardsson et al 

iReddv et al. 

'^1 

ol 

Casagrande et al. 


[1993b; lRocha|pint^ e^ ni2QQ0l); 
i2Qllb ; [Bensb^Feltzing. & OevI 


blue full triangles, magenta full squares, grey small dots, green stars, and 
purple crosses, respectively. The large yellow dots are the solar neighbour¬ 
hood SFR at the present time in a) and the solar neighbourhood metal- 
licity at the time when th e Sun formed, 4. 5 Gyr a go, in b). c) [(r/Fe] as 
a function of [Fe/H] from ICasagrande et alj <201 ih (cyan dots). The large 
yellow dot represents the solar abundances at coordinates (0,0). d) The 
MDF with data from Chang. Hou. & fJ i2QQ(l) ; ICasagrande et alj j2Ql ih 
and iKordopatis et" i] boiSh tRAVE survey) as magenta squares, green tri¬ 
angles and black dots, respectively. In all panels the red dots with error bars 
are the binned results using all the noted datasets. 


occurred at a time 3-5 Gyr after t - 0, in agreement w ith more re¬ 
cent works of icignoni et alj j2Q06h an d[R owe 3 (l2013h . These data 
are binned in 1 Gyr time-steps for the analysis which follows. The 
results, given in Table I All have then been normalised to the most 
recent values of the SFR in the solar region, corresponding to the fi¬ 
nal point at 13.2 Gyr. The value for the present-day SFR for the en¬ 
tire MWG is estimated to be in the range [0.8-13] MQyr“* jRanaL 
Misiriotis et al. j2006h give a value of 2.7MQyr“*, while 


Chomiuk & Povichl j201 ih find 1.9MQyr *. Taking into account 


the ratio of areas of the MWG disk within the optical radius and 
that of the solar region, we obtain a value of 'Tq ~ 0.266 Mq yr“' 
for the region located at a galactocentric distance R = 8 kpc, in ex¬ 
cellent agreement with the used value in ICalura et alj ( 1201 dl . This 
value is the large yellow dot shown in Fig. lAlb l. In this figure, 
we see the time evolution of the SFR once normalised to recover 
this value Tq ~ 0.266 at 13.2 Gyr. The values of T are given in 
Table IaTI 

In panel b) of Fig. lAll we show the age-metallicity rela¬ 
tion obtained with data from the literature as labelled, binned for 
each Gyr, as in pane l a). Data from recent surveys suc h as RAVE 
dSoeche et al.Ll20f3) or APOGEE l l Anders et all l20l3l fall in the 
same region of the plane in panel b) when only the solar regioiFI 
data are selected. Our binned results are shown in Table IaTI They 
have been normalised to obtain a value [Ee/H]=-t0 in R = 8 kpc 
at a time t = 8.5 Gyr, when the Sun was bom, implying a shift 
of -1-0. Idex compared with the data of Fig. lAlb l. In both cases. 


The solar region is defined as a 1 kpc annulus centered on a galactocen¬ 
tric radius of 8 kpc, with a thickness of 200-500 pc. 


we have added to the dispersion obtained from the binning pro¬ 
cess, a systematic error (representing observational uncertainties) 
of 0.05 dex and O.lOdex in columns 3 and 5, respectively. In panel 
c) we show the values of the a-ele ment abundances compar ed with 
those of iron, with data taken from ICasagrande et al.l ( 1201 ill . From 
the latter, we have selected those stars located between 7.5 and 
9.5 kpc that lie within 0.5 kpc pf the mid-plane of the disc, for 
studying the evolution of the solar neighbourhood. These values are 
binned and shown in Table IA2I In panel d) we show the metallic- 
ity distribution function (MDF) . Given the similarity of the three 
datasets, we have binned and normalised the result to unity, and 
listed them in Table lA2] 


A2 The MWG disk: surface densities 


The radial gas distributions for both molecular and diffuse phases 
are well known. Since our model calculates separately both compo¬ 
nents, we also use these observations to fit our models. We use data 
from the literature, shown in Fig. lA2l for diffuse HI and H2. By bin¬ 
ning both sets, we obtain the results given in Table lA^ and shown 
in panels a) and b). We see clearly a maximum around 10 kpc for 
HI while H2 shows an exponential shape from 4 kpc to the outer 
disk. It also shows the well known molecular hole inside ~3 kpc. 

In panels c) of the same Fig. lA2l we also show the stellar sur¬ 
face density profile, including estimates from different authors as 
labelled. The most recent estimates for the solar stellar surface den- 


siiy give values Deiween ana 04 pc ” aiyuiiKen & i.ji 

19891 I1991I; IVallenari. Rertelli. & Schmidtobreickl 

mora 

20Qd; 

Siebert. Bienavme. & SoubiranL l2003l: iKhooerskov & TvurinaL 

2OO3I; iHolmberg & Flvnnl 20041; [Bienavme et all 

20061; 

Flvnn et all I2OO6I: IWeber & de Boeil l20 C 

; iMcMillanl 

201 ll; 

Moni Bidin et all 20121; ISurch & Cowsik 

2OI3I; IZhang et al.l 


20131) . These values depend on the scalelength 


for the disk Rj, which is in the range [2.15^] kpc. We show these 
data in panel c) of Fig. lA2l with our results after binning (red dots). 
In panel d), we show the SFR normalised to the solar value . The 
binned results for each kpc are also shown as red points, as in 
panels a), b), and c). 

In Table IA3l we present the resulting binned-averaged values 
of diffuse and molecular gas surface densities, and their associated 
errors, (columns 2 to 5), in Mq pc“^, for each radius given in col¬ 
umn 1. The stellar surface density profile is given, in logarithmic 
scale, with its associated error, in columns 6 and 7. In columns 8 
and 9 we show the SFR surface density in Mq pc“^ Gyr“*. 


A3 Disk elemental abundances for C, N, and O 

C, N and O abundances are the most important constraints for our 
models. Since N comes mostly from intermediate mass stars, O 
from the massive ones, and C from both, a fine-tuning of the stel¬ 
lar yields and IMF is necessary to reproduce simultaneously the 
three elements. We hope that any of the different combinations of 
yields from low and intermediate mass stars, with those from mas¬ 
sive ones, with different IMFs, would give the right CNO elemen¬ 
tal abundances. We show in Fig. IA3l the three radial distributions 
for C, N, and O in panels a), b) and c), respectively. Data from 
different studies are plotted with different symbols, as listed in Ta¬ 
ble |A4l while the red large dots are again our binned results (as 
12-1- log {XjH)) shown in Table I A3] 
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Table A3. Radial binned distributions obtained from observational data 


R 

(kpc) 

^Hl 

error 

error 

M© 

logs. 

error 

iog'E.SFR 
M© pc“2 

error 

Gyr* 

C/H 

AC/H 

N/H 

AN/H 

0/H 

AO/H 

0 

9.41 

1.00 

0.30 

0.50 



-0.370 

0.15 

(8.66) 

(0.30) 

8.39 

0.31 

9.02 

0.40 

1 

3.97 

1.88 

3.82 

4.90 



0.603 

0.52 

(8.64) 

(0.30) 

(8.24) 

(0.30) 

(8.86) 

0.30 

2 

2.37 

2.04 

5.18 

5.39 



0.706 

0.47 

(8.55) 

(0.30) 

(8.20) 

(0.30) 

(8.74) 

0.30 

3 

2.39 

2.12 

3.48 

2.24 

2.43 

0.01 

0.983 

0.59 

8.35 

0.67 

7.84 

0.59 

8.62 

0.45 

4 

3.86 

2.35 

5.69 

3.35 

2.50 

0.13 

1.163 

0.35 

8.63 

0.15 

8.16 

0.29 

8.82 

0.45 

5 

5.06 

2.14 

8.28 

2.87 

2.40 

0.07 

1.185 

0.24 

8.48 

0.27 

8.02 

0.66 

8.83 

0.33 

6 

5.04 

2.06 

8.47 

1.67 

2.25 

0.07 

1.181 

0.27 

8.40 

0.32 

7.83 

0.36 

8.77 

0.56 

7 

5.44 

1.58 

4.59 

1.72 

2.09 

0.08 

0.963 

0.26 

8.34 

0.25 

7.87 

0.41 

8.69 

0.30 

8 

5.69 

2.38 

3.15 

1.42 

1.95 

0.09 

0.123 

0.29 

8.28 

0.20 

7.77 

0.30 

8.56 

0.30 

9 

7.69 

2.13 

2.44 

0.80 

1.79 

0.10 

0.594 

0.25 

8.26 

0.22 

7.76 

0.41 

8.60 

0.35 

10 

6.52 

2.18 

1.96 

1.18 

1.69 

0.14 

0.510 

0.43 

8.11 

0.27 

7.59 

0.37 

8.45 

0.36 

11 

6.16 

2.06 

1.24 

0.80 

1.51 

0.03 

0.403 

0.39 

8.01 

0.29 

7.60 

0.36 

8.41 

0.34 

12 

5.63 

2.17 

0.99 

0.75 

1.38 

0.01 

0.006 

0.65 

8.00 

0.38 

7.53 

0.31 

8.44 

0.32 

13 

4.83 

2.86 

0.57 

0.71 

1.25 

0.03 

0.183 

0.48 

7.76 

0.20 

7.35 

0.33 

8.44 

0.38 

14 

3.65 

2.80 

0.82 

0.94 

1.09 

0.01 

-0.260 

0.57 

7.93 

0.17 

7.45 

0.34 

8.42 

0.43 

15 

2.96 

2.69 

1.09 

1.84 

0.94 

0.01 

-0.132 

0.59 

7.60 

0.16 

7.36 

0.47 

8.14 

0.43 

16 

2.42 

2.19 

0.20 

0.07 

0.80 

0.01 

-0.520 

0.15 



7.60 

0.54 

8.14 

0.39 

17 

2.15 

2.17 

0.13 

0.05 



-0.680 

0.18 



6.98 

0.62 

8.19 

0.35 

18 

1.61 

1.77 

0.08 

0.03 



-0.890 

0.15 





7.96 

0.50 

19 

1.18 

1.66 

0.03 

0.01 



-1.370 

0.15 







20 

1.10 

1.60 
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Figure A2. Radial distributions of surface densities in the MWG for: a) 
diffuse gas, Z/// (in Mppc " ^ units) with data from lOlling & Merrifielc 
2001 ): [Wolfire et ^ i2003 ():|Nakanishi & Sofuj ^200^ : iKalberla 


200 ^ ) and IPineda et alJ i2013h . as cyan asterisks, orange stars, blue full 


triangles, green full dots, and magenta full squares, respectively; b) molec- 


ulai'gas, Ljy, (inM(^pc ^ u nits) with data fromlwilliams & McKee 
iNakanishi & Sofu'3i2006h : |Pineda et alJ i2013h and lUrquh^^^l 


1997): 


'20141) 


as green full triangles, blue full dots, magenta full squares, and cyan 


Talbotlh980l);lRand Jl991 

);IVallenari. Bertelli. & SchmidtobreickI i2000l): 

Bow & Rixl j2013l);ISofue 

(20131); d) the SFR surface density. 

normalized to the Solar 

value = 0.266 Mrr>vr“*, estimated from 

Misiriotis et alj j2006h;IChomiuk & Povichl J201 ih in R = 8kpc, with data 

fromiLacev & Falll il983h;lwilliams & McKed ^997|). as ureen full trian- 


pemovae and Hii regions, shown as yellow stars, orange open dots, an d 
magenta full squares. We have also used those from IUrquhart et ^ j2Q14h , 
represented by cyan asterisks. These last two panels are in logarithmic scale. 
The binned results ai‘e the large red dots with error bars in all of them. 



Radius (kpc) 


Figure A3. Radial distributions of abundances (as 12 + log {XIH)) for a) C, 
b) N and c) O. Data are taken from the works noted in Table lA4l where the 
symbol used for each is also given. In all panels the red full dots with error 
bars are the binned results obtained in this work and given in Table lA3l 
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Table A4. List of data sources employed in Fig. |A3l 


Author 

c 

N 

o 

Symbol 

Peimbe^^^^) 

- 

X 

X 

black * 

Shaver e^al_^983) 

- 

X 

X 

orange o 

Fich & siiig£ ywn 

- 

X 

X 

yellow X 

Vilchez & Estebai^^96) 

- 

X 

X 

green ★ 

Afflerbach. Churchwell. & Werner (1997) 

- 

X 

X 

blue ■ 

Esteban. Peimbert. & Torres-Peimbert (1999) 

X 

X 

X 

green o 

Reddy et al. f2003) 

X 

X 

X 

blue A 

Daflon & Cunha (2004) 

- 

X 

- 

magenta o 

Esteban et al. (2005) 

X 

X 

- 

brown □ 

Gavilan. Molla. & Buell (2006) 

X 

X 

X 

black • 

Rudoloh et al. (2006) 

- 

- 

X 

magenta o 

Henrv etaL_^20W) 

- 

- 

X 

blue • 

Balser et^j_^20^) 

- 

- 

X 

green A 

Luck e^^^2^n ) 

X 

X 

X 

cyan ■ 

Esteban et al. (2013) 

X 

- 

- 

green A 
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